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Abstract 



Assuming the Hubble parameter is a continuous and differentiable 
function of comoving time, we investigate necessary conditions for 
quintessence to phantom phase transition in quintom model. For power- 
, law and exponential potential examples, we study the behavior of dy- 

' namical dark energy fields and Hubble parameter near the transition 

time, and show that the phantom-divide-line lo = — 1 is crossed in these 
' models. 

o 

1 Introduction 

o 



Astronomical data show that the expansion of our universe is accelerated at 
the present epoch pQ. Assuming that the universe is filled with perfect fluids, 
the equation of state parameter uj = p/p must satisfy u < — |, indicating a 
negative pressure. Many theories have been proposed to study the origin of 
this negative pressure or repulsive gravitational behavior. 
■ One of these theories introduces a smooth energy component with neg- 

ative pressure dubbed "dark energy". A candidate for dark energy is the 
cosmological constant ;2, a : a constant quantum vacuum energy density which 
fills the space homogeneously, corresponding to a fluid with a constant equa- 
tion of state parameter uj = — 1. Because of conceptual problems associated 
with the cosmological constant, such as fine-tuning and coincidence prob- 
lems 0, alternative theories have been proposed where in a class of them, 
some dynamical scalar fields, with suitably chosen potentials, have been in- 
troduced to make the vacuum energy vary with time [2]. For example to 
describe u > — 1 phase or quintessence regime, a normal scalar field 4>, 
known as quintessence scalar field, can be used @]. For performing actual 
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calculation, it is useful to assume that the quintessence field is slowly rolling 
in a potential. This approximation (slow-roll), leads to inflationary expan- 
sion of the universe • u) < — 1 phase (or phantom phase) can be related to 
the presence of a scalar field a, with a wrong sign kinetic term dubbed as 
phantom scalar field Depending on the form of the potential, different 
solutions such as asymptotic de Sitter, big rip, etc. may be obtained [Zj. The 
effects of gravitational back-reactions can also counteract that of phantom 
energy and can become large enough to end the phantom dominated phase 
before the big rip [H]. 

One of the important issues of the scalar models of dark energy is their 
stability behaviors in both the classical and quantum mechanical levels. In 
classical level, one of the main methods of stability studies is achieved by 
studying the stability of the late-time attractors of the theory by a phase 
space analysis. This can be done by determining the eigenvalues of the de- 
terminant obtained by the set of autonomous equations of motion when 
perturbed about their critical points. Variety of models have been studied 
by this method, which some of them can be found in [§| and references 
therein. 

The second kind of instabilities, which is more important in phantom 
models, are those based on the quantum fluctuations. Because the phantom 
fields have negative kinetic energy, it is possible that a phantom particle 
decays into arbitrary number of phantoms and ordinary particles, such as 
gravitons. It can be shown that the decay rates of these interactions are 
infinite which indicates that the phantom models are dramatically unstable. 
But if we think of these models as the low-energy effective theories, with 
the fundamental fields having positive kinetic energy, then we should use 
a momentum cutoff A in calculating the decay rates. In this way it can be 
shown that for A ~ M p \, the lifetimes can become larger than the age of the 
universe when one chooses suitable phantom-gravity interaction potentials, 
and this removes the quantum instability of these kinds of phantom models. 
See ^U] for a specific example. 

Some present data seems to favor an evolving dark energy with uj less 
than —1 at present epoch from uj > —1 in the near past It is 

therefore instructive to construct physical models in which uj can cross uj = 
— 1 line (dubbed as phantom-divide-line). Neither the cosmological constant 
nor the dynamical scalar fields, like quintessence or phantom, can explain 
the uj = —1 crossing, hence some other theories have been introduced to 
describe this transition JH]- One of these theories is the quintom model, 
which based on hybrid models ^T5\, describes the transition from uj > — 1 
to w < —1 regime by assuming that the cosmological fluid, besides the 
ordinary matter and radiation, is consisted of a quintessence and a phantom 
scalar field A brief introduction of quintom model can be found in 
section two. This model can lead to quintessence domination, i.e. u > — 1, 
at early time and phantom domination, i.e. uj < —1, at late time. In |16j . a 
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phase-space analysis of the evolution for a spatially flat Friedman- Robertson- 
Walker (FRW) universe containing a barotropic fluid and phantom-scalar 
fields with exponential potentials has been presented. It has been shown that 
the phantom-dominated scaling solution is the stable late-time attractor. In 
|17j the same calculation has been done by introducing an interaction term 
between phantom and quintessence fields. Recently, the hessence model, as 
a new view of quintom dark energy, in which the dark energy is described 
by a single non-canonical complex scalar field, rather than two independent 
real scalar fields, has been introduced |18|. The evolution of u; in this model 
has been studied in JH] and [201, y i a phase-space analysis. Although this 
model allows u to cross —1, but it avoids the late time singularity or the 
"big rip" . 

In this paper we study the transition from quintessence to phantom 
phase in the quintom model. Instead of considering late time behavior of 
the model we try to investigate the behavior of cosmological parameters 
and scalar fields in a neighborhood of transition time. In section one, by as- 
suming that the Hubble parameter is continuous and differentiable function 
at transition time, we seek the necessary conditions that must be satis- 
fied in order that the phase transition occurs. In section two, instead of a 
phase-space analysis, we try to obtain solution of Einstein and Friedmann 
equations in the vicinity of transition time for exponential and power law 
potentials for slowly varying fields. We show that the dynamical equations 
are consistent with the conditions needed for quintessence to phantom phase 
transition. It is worth noting that these two specific examples have stable 
late time attractor solutions ^§1 El • 

Through this paper we use h = c = G = 1 units. 



2 Phase transition in quintom model: general re- 
sults 

We consider a spatially flat FRW space-time in comoving coordinates (t, x, y, z) 

ds 2 = -dt 2 + a 2 (t)(dx 2 + dy 2 + dz 2 ), (1) 

where a(t) is the scale factor. We assume that the universe is filled with two 
kinds of fluid: (dark) matter and quintom dark energy. The equation of state 
of matter is 

Pm = (7m - l)Pm, 1 < 7 m < 2, (2) 

where p m , and P m are the matter density and pressure, respectively, and 
7 m = 1 + uj m - The evolution equation of p m is 

p m + 3Hj m p m = 0. (3) 
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H(t) = a(t)/a is the Hubble parameter and "dot" denotes time derivative. 
The quintom dark energy consists of a negative kinetic energy scalar field a 
and a normal scalar field <fi, described by the Lagrangian density [TT] . [T6] : 



Ld = \d^<t> - ^ad^a + V{<t>, a), ( ! ) 



where V(cp, a) is the quintom potential. Restricting ourselves to homoge- 
neous fields, the energy density pr> and pressure Pd of the homogenous 
quintom dark energy is then 

PD = \ft ~\o 2 + V(<j>,a), 

Pd = \fi-\d 2 -V^o). (5) 
The evolution equations of the fields are 

ocp 

a + 3Ha-^^l = 0. (6) 
oa 

Using Einstein equation, one can show that the Hubble parameter satisfies 
the Friedmann equations 

H = — {PD+Pm) 

= —[^-a 2 + 2V(^a) + 2p m ], (7) 

and 

H = -4vr(p D + Pm + P D + P m ) 

= -4tt(<P 2 - a 2 + 7 m /j m ). (8) 

Note that eqs.((SJ), (0 and, (jHJ are not independent. The equation of state 
parameter uj is defined through 

Pp + Pm = 4> 2 -a 2 - 2V(</>, a) + 2( 7m - l)p m 
PD + Pm <j) 2 -& 2 + 2V(<t>, a) + 2p m 

In terms of Hubble parameter we have 

2 H 

" = -l"3^- (10) 

For an accelerating universe, d(t) > 0, we have H(t) + H 2 (t) > and 
u> < — |. Through this paper we will restrict ourselves to H(t) > 0. In the 
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quintessence phase , H < and therefore to > — 1. In the phantom phase 
we have H > so w < —1. In contrast to cosmological models with only 
one scalar field, transition from quintessence to phantom era, in principle, 
is possible in quintom model. If H(t) has a local minimum at t = to, i.e. 
H(t ) = 0, then at t < t , H(t) < and u) > -1 and at t > t , H(t) > 
and u < — 1. This behavior is in agreement with present data of equation of 
state parameter u. So the transition from quintessence to phantom phase, 
or crossing the phantom-divide-line, can be studied by investigating the 
behavior of the Hubble parameter near t = to point. 

Let H(t) be differentiable on an open set containing to, then it can be 
expanded at t = to 

H(t) = Y, T^-*°) n ' H(t ) = 0, (11) 

where H^ n \to) is the n-th derivative of H(t) at t = to- If a > 2 is the order 
of the first non-zero derivative of H{t) at t = to, then 

H(t) -ho + hit- to) a + 0((t- t ) a+1 ) . (12) 

where ho = H(to) and hi = -^H^ a \to)- A transition from quintessence 
to phantom phase occurs when a is an even integer and hi > 0. hi > 
follows from the fact that H(to) must be a local minimum of H(t). If we 
considered the transition from phantom to quintessence universe, we should 
take hi < 0, instead. Using eqs. 0, ((HJ), one finds 

H " = f V " ^' (13) 
where V(t) = V((j), a) + (1 — ^)p m - By expanding both sides of eq. ((T8|) near 
t = to, and noting H^ <a \to) = 0, one finds 

V^-^io) = 0, (14) 



and 



hi = f V , (15) 

h = —Vi, (16) 

a 

in which Vo = V(to) and Vi = t^zw V^ a ~ 1 \to). For matter density p m , one 
P&T 1} (*o) = (SHitoh^Pmito), (17) 



finds, using eq.© 



so 

kl = f[ (v^Wo)'^)) + (1 - ^ i )(-37 m / i o) Q -V m (to) j • n > 2. 



(18) 
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Assuming p m (t) > and noting j m < 2, the second term in the right-hand- 
side of ea.([18|) is negative for even a's. So to have hi > 0, it is necessary 

W a - 1 )(0(t o ),cr(t o ))>O. (19) 

It is also interesting to note that ea. (|15|) implies 

V(to) > 0. (20) 

But V(t) > is a necessary condition for acceleration of the universe at any 
time t. This can be verified by considering that when uj < — |, eq.@ yields 

ft -d 2 < + (l - ^7m) Pm- (21) 

Also for H(t) / 0, eq.(J7J) results 

4> 2 -& 2 > -2V{^a)-2 Pm . (22) 

V(t) > follows from ens. lf2Tj) and (1221 . 

Let us for a moment relax our even-a condition and look at the models 
with odd a. Assume that at some time t we have H(t = t) = 0. If the system 
is in quintessence phase in t < t times, i.e. H(t < t) < 0, expanding H(t) 
near t = to results 

H(t) = h + h!(t - t) 2n+1 + 0((t - i) 2n+2 ), n > 1. (23) 

We assume hi < to guarantee the existence of quintessence phase in t < t 
times (for hi > 0, the system describes a phantom universe). In this case t 
is an inflection point and H{t) < for both t > t and t < t, and therefore no 
transition to phantom phase happens. Assuming the potential has a lower 
positive bound V((p, a) > vo > 0, and using the fact that H{t) < for all 
t's, eqs.® and (jHJ) result 

H 2 (t) > 

> 

Since H{t) is a decreasing function, the above relation shows that H 2 (t) 
achieves its minimum at infinity ( provided we assume that the system re- 
mains in quintessence phase for all t) . Thence lim^oo H 2 (t) = (87r/3)^o, 
while lim^oo V(cp, a) = vq and lim^oo p m = 0, and the system tends to a 
de Sitter space time at late time, i.e lim^oo H{t) = |21j . 

Now let us come back to even-a case and study the behavior of the 
equation of state parameter u near t = to, where H{to) = 0. For non-zero 



47T 

8vr 



((2--/ m ) Pm + 2V(cP,a)) 



-Vq. 



(24) 
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differentiable H{t), u can be Taylor expanded as (to see a discussion about 
Taylor expansion of oj in terms of scale parameter see |22| ) 

^) = E^T M (^^or. (25) 



n=0 



TV. 



From eq.(|10[). it is found that a; (to) = — 1 and that the first non-zero deriva- 
tive of u)(t) in its Taylor expansion is 

, ,(a-l)/ + \ _ 2 g(a) (^o) / 9R n 

U (to) -"3^oT' (26) 

where as before we have denoted the order of the first non zero derivative 
of H{t) at t = to by a. Therefore 

w(t) = -1 - f - ^o)"- 1 + 0((t - t r), (27) 

or in terms of V, 

u(t) = -l-^(t- to)"- 1 + 0((t - to) a ). (28) 

The above equation shows that (i) to cross the phantom-divide- line w = — 1 
a must be an even integer, (ii) to go from quintessence phase, to > —1, to 
phantom phase, u> < —1, we must have h\ > 0, or equivalently Vi > 0. Note 
that if initially the universe is in phantom phase, the condition (ii) for phase 
transition to quintessence era becomes h\ < 0. These results are consistent 
with those we have obtained from studying the behavior of H (t) near t = to- 
In brief, to have an accelerating universe, a specific combination of the 
quintom potential and matter density, denoted by V(i), must be positive. 
Also a transition from quintessence to phantom phase occurs provided the 
order of the first non- vanishing derivative of H{t) at t = to, in which 
H(to) = 0, is even and the parameter hi obtained in eq.(|18|) is positive. 
We will illustrate these general results and discuss the phase transition via 
two important examples: exponential and power like quintom potentials, in 
the next section. 



3 Examples 

3.1 Exponential potential 

As a first example we consider the quintom model with potential 



v = v d + v a 



vie 



Ai0 



+ v 2 e Ma ,vi > 0, v 2 > 0. (29) 
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To study the possible occurrence of transition from w > -1 to w < -1 for 
this potential, we try to obtain the solutions of eqs.©, 0, an d ©j when 
the Hubble parameter behaves as 



H(t) = h + M Q + 0(t a+i ), a>2,hi^0, 



(30) 



near t = 0, in which H(0) = 0. Eq.lJHJ) cannot be solved exactly with the 
potential (|29|). but for slowly varying fields, where <j> « Hep and a « Ha, 
these equations become 



3(/to + h x t a )^{t) = -XiVie Xl<t>{t) 
3(h + h x t a )a{t) = \ 2 v 2 e x ^ t \ 



(31) 



with solutions 
a(t) 



3/ioa 



Ai | UlA 2 td>(^Ai,i,I) +Cl/loa 



(32) 



3/iqo 



W2A 



ho 



1, +c 2 /i a 



in which $(2, a, 6) is the Lerchphi function. The constants c\ and c 2 , in terms 
of the initial conditions at t = 0, are defined through 



1 



■In 



v 1 Xi ci 



(33) 



(7(0) = r ln 



V2A2 C 2 



Hence c\ > and C2 < 0. Near i = 0, up to order O(t), we have 

0(0) 1 

3/100(0) 



^cf Ai ' 
3 

Aici ' 



(34) 



and 



*(0) 
3h &(0) 



1 

3 



A2C2 



(35) 



Therefore the slowly varying condition is equivalent to 



\cihl\ » 1 , |c 2 /iol > 1- 



(36) 
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These inequalities can be expressed in terms of the potentials as following: 
li(^) 2 (0)l « h 2 and |^(^) 2 (0)| « ftg, respectively. 
p m (t) satisfies the equation 

Pm{t) + 3 7m (/i + hxt a )p m {t) = 0. (37) 

The solution of this equation is 

Pmit) = c 3 e- 3 ^ + ^ +1 \ (38) 

where c 3 = p m (t = 0). 

In terms of dimensionless variables C\ = 1/(ci/iq), C2 = l/^/t 2 ,), C3 = 
C3//10) r = ^oi; an d H\ = /ii//iq +1 , eq.® becomes 



H 2 (t) _ 4vr 



a 



r*(-r«ffi, 1, ±) + #1 (l + tfxr^A 2 



(39) 



cr 



+ 



6a 



^(-^,1,1)+ al 2 (1+jffirQ) 2 A 2 ' A 2^ ( _ T a ifl)1; I ) + ^ 



6a 



\ 2 
A 2 



^(-^,1,1) + ^ 



+ 2C 3 e~' 37m ^^ 



and eq. (JHJ) reduces to 
1 dH{r) 



ho dr 



4ir< 



a 



T^-T^X, + (1 + //IT") 2 A? 



(40) 



r*(-T«fli, 1, ±) + ^1 (1 + tf lT <*) 2 A 2 



By expanding the right hand sides of eqs. (|39|) and (f4T)|) near r = and using 
eq.fjSOJ) for their left hand sides, one finds: 



1 + 2H lT a + 0(t 



Air 



Ci(Ci + 6) C 2 (C 2 + 6) | 



A 2 



A 



-8vr 



47T 



cm + 1) cKf+i) 

- - - r 7mW 



A? 



A| 



T 



3 [ 3 7 2 t C 3 + 2^ - 2^f ] +21^-^1^ 



+3 



A 2 A| 
r 2 + 0(r 3 ), 



(41) 
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and 




/ c 2 
-4tt I 7m C 3 + 



2Cf 



+ 47r|3 7 ^C 3 + ^-^lT + 0(T 2 ), (42) 



respectively. To have a consistent set of equations for a > 2, one finds the 
equalities 

4vr /Ci(Ci + 6) C 2 (C 2 + 6) 



and 



3 V Af A| 



C 2(f + 1) Cf (f + 1) 



Af X 2 



+ 2C 3 = 1, (43) 



+ 7mC 3 = 0, (44) 



from the first two terms of eq. (|41[) and 

lmC 3 + ^ " = 0, (45) 



Cf_C? 



from the first term of eq.(|42|). respectively. In slow-roll approximation, which 
we are working in, the parameters satisfy ea. ()36|) . which imply \C\\ « 1 
and I C2 1 << 1. In this approximation, the eqs. (|44j) and ()45jl are equivalent. 

The important observation is that since C\ and C3 are positive and C2 
is negative real numbers, the coefficient of r in eq. (|42|) is a non-zero positive 
number, which results the following for exponential potential: 

C^exp. — 2, 

(Mcxp. > 0. (46) 

In this way it is proved that in the quintom model with exponential potential, 
the system has a phase transition from quintessence to phantom phase, if any 
other remaining relations, obtained from eqs. (|41|) and (|42|) . will be satisfied 
consistently up to the lowest order. 

To check the remaining relations, we first note that since a = 2, ea. (J42j) 
results another equation, as follows, 



2,(37^3 + ^-^]. (17) 



But eqs.jSI) and gHJ imply: Cf/Af - C|/A| « \C 2 /Xj - Cf/A|| ~ C 3 , so 
eq.(|47|l reduces to: 

Ht ~ 67r 7 ^C 3 . (48) 

This shows that the rate of the phase transition depends on the parameter 
of state j m and the dark matter density at transition time. 
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In terms of dark energy and matter densities, po, Pm, eq . (|43j) can be 
written as 

— + — = 1, (49) 

Pc Pc 

where p c = 3/ig/8vr, p m (0)/p c = (8^/3)C 3 and p D (0)/p c = (4tt/3) [d (d + 
6) / Af — C 2 (C 2 +6)/ X^]- The dark energy density is composed of two parts: the 
density of kinetic dark energy, (47r /3)(C 2 / X 2 — C 2 /X 2 )p c , which is negative 
according to ea.(|43|). and the density of potential dark energy, (47r/3)(6Ci/Af - 
QC 2 /X 2 2 )p c . Now since Ci/Xj - C 2 /X 2 2 » C 2 /Xj - Cf /Al ~ C 3 , it is obvious 
that the main part of the energy density at the transition time is coming 
from the potential of dark energy. 

The remaining equation that can be extracted from ea.(|41|) is 



Hi 




+ 2 



A?, 



Hi 



+ 0(C 3 ) 



(50) 



which can be verified using the fact that Cf/X 2 — Cf /A| << Ci/X 2 -C 2 /X 2 ~ 
Ptotai/Pc = 1- This completes our desired checking of all relations in the 
lowest order. 

Expanding V(t)/hl = (V + (1 - ^)p m )Ao at t = results 



hi 



(2-7, 



]C3 + 3Ci 



+: 



A? 



3C 2 
A| 
2C| 
A! 



+ 



2Cf 
A? 



r + 0(r 2 



(51) 



where we have used eas.(|43)l. and (|15]l in the second equality. The 

above result is in complete agreement with our previous assertions eqs. (|15j) 
and HI}). 

To test, in another way, the validity of our approximations let us choose 
Ci = 10 -4 , C 2 = —1.1 x 10 -4 ( since it must be small quantities ), "y m = 1, 
and Ai = X 2 ( for simplicity ), at t = 0. Eqs.(jJ3), (0 and (JEJ) then result 
C 3 = 4.0 x 10- 7 , Ai = A 2 = 0.073, and Hi = 7.5 x 10~ 6 . Fig.© shows the 
plot of (1 + Hit 2 ) 2 and the same plot obtained from ea. ()39|) . Fig.© shows 
the comparison of to in these two approaches. 

At the end, it is worth noting that all the above approximations are valid 
until r << 1, which is equivalent to t « . But this is reasonable period 
of time since 1 is of order of the age of our universe. 
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1 .00002 



Figure 1: H 2 /h^ as a function of r, obtained from ea. (|39|) (points), and 
(1 + tfiT 2 ) 2 (line). 
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3.2 Power-law potential 

As the second example consider the potential 



V 



Vs + Vv 



= v 1( j> bl + v 2 a b \ 
Using the approximation (|30j) . the solutions of the eq. 



(52) 



are 



0(t) 



a(t) 



'v 2 b 2 (-b 2 + 2)i$(- 



+ ci 



1,4: 



i 

2-b-i 



3/iqQ; 



+ C 2 



(53) 



where a = 2 ~ bl (O) and c 2 = cr 2 ~ b2 (0). We assume 0(0) > and cr(0) > 0. 
Using these solutions one can show that neglecting second derivatives in 
eq.©, is allowed when 



|6 1 (6 1 -1)^|^- 2 (0)<1, 
|6 2 (6 2 -l)§k fe - 2 (0)«l. 

i d 2 V rh 



(54) 



In term of potentials, these inequalities are |^£|(0) <C h\ and |^t"|(0) <C 
h\. For b\ = b 2 = 2; v\ = m 2 /2 and v 2 = m 2 /2, ea.l|53|) reduces to 



0(t) 



0(O)e" 



3/lQQ 



<r(t) = a(0)e" 



3/lQ a 



(55) 



In this case, the approximation (|54[) is reduced to small mass limits for 
phantom and quintessence fields: <C /io and mo- -C ho. 

p m is given by eq.([38|). By putting eqs. (|38|) and (|53]> into eqs.® and 
©, and using dimensionless variables Vi = Vi/ho, V 2 = v 2 jh\, C3 = c^/h^, 
t = tho, and Hi = hi/h$ +1 , near r = 0, we obtain 



4-7T 

1 + 2H X T a + 0(r Q+1 ) = — 



9 



2F 2 a b2 (0) 
2 



47T 

+ T 



26? (61 - l)Xf0 36l - 4 (O) 2&3(6 2 _ l)U 2 3a 3fe2 " 4 (0) 



27 



27 



67mC 3 
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+- 



4tt 



27 



+ 



^2( 6l -l) (0)6 2y2 a 2{b 2 -l) (0)6 2y s 



+ 



27 



#1 + 



6f (36? -76i +4)y 1 V 4fel " 6 (0) 



81 



6f (36| - 76 2 + 4)F 2 V 4fe2 - 6 (0) 



81 



r 2 + 0(r 3 



(56) 



and 



a.H\T a + 0(t c 



-47T 



9 



9 



+ 7m C 3 



-47T 

+0(r 



26 3 (6! - l)0»i- 4 (O)l? 26|(6 2 - l)a 3 ^- 4 (0)y 2 3 ' 

-37m C 3 



27 



27 



(57) 



respectively. Therefore if the assumption (|30|) is true, we must have 



^2(61-1) ( )&2^2 a 2(6 2 -l) (0)6 2 F2 



2 +7mC 3 = 0, 



(58) 



7mC 3 + - 



(Q) A + l 6l(6l _ 1)^1-2(0) 



V%bla 2 ^- l \<d) (l - h 2 (b 2 - l)V 2 ^- 2 (0) 



0.. 



(59) 



and 
4tt 

T 



„ )W ^-(OW^ + 2Ca + 2Fi/l(0) + ^ ' 



9 



9 



1. 



(60) 



For slowly varying fields (eq.(|54|)). eq. (|59j) reduces to eq. (|58j) . The coefficient 
of r in ea, (|5Tl) is 



47T 



r26? (bi - l)^ 361 " 4 ^)^ 3 2b 3 2 (b 2 - 1)ct 362 - 4 (0)T/ 2 3 



27 



+ 



27 



(61) 



Comparing eqs.® and (JHHJ shows ^ 3bl - 4 (0)V 1 3 +cj 362 ~ 4 (0)y 2 3 << ^^^^O)^ 2 - 
cj 262 ^ 2 (0)F 2 2 ~ C 3 , which reduces eq.JHU) to 



#i ~ 67T7 m C 3 . 



(62) 
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This is a non-zero positive number. So for quintom model with power-law 
potential, we show 



empower— law — 2, 
(frl)power-law > °> ( 63 ) 

which proves the existence of quintessence to phantom phase transition in 
this model. Of course the remaining relations must be also checked. Note that 
the value of Hi at transition point does not depend on the potential, compare 
eas.(|48|) and a fact that we guess is true for any other potential. 

To check the consistency of the remaining equations, first we note that by 
the same arguments used in the preceding example, it can be shown that the 
main contribution in ea. (|6(Jj) . which can be written as ptotai/Pc = 1, comes 
from the quintom potential, i.e. Vi(j) bl (0) + V2«r b2 (0). Now the coefficient of 
t 2 in the right-hand-side of ea, ()56|) . using ea. (|oTj) . can be written as 



2H 




^(fti-l) (Q) 6 2 F 2 ^2(62-1) (Q)fr2y 2 2 \ 



9 9 



(64) 



which is equal to 2H 1} using the fact that -V^cf 2 ^ 1 '^ (0)+V 2 2 a 2( - b2 ~ 1) (0) << 
Fi0 bl (O) + V 2 a b2 (0) ~ ptotai/Pc = 1- Finally it can be shown that eq.(JHU) is 
also verified for this potential. 

As an illustration of the transition behavior of quintom model with 
potential (JS2J), see Fig(3), which shows the plot of H 2 /Kq obtained from 
eq.0, for b x = b 2 = 2, -y m = 1, vi = m^/2, v 2 = ml/2, with m<j>/h = 
ma /ho = 0.01 ( since they must be small quantities ) and <7 2 (0) = 2<p 2 (0). 
Eqs. (|SSJ, dSOJ) and (ggj then result ^ 2 (0) = 2500.0/tt, C 3 = 0.28 x 10- 5 /tt, 
and Hi = 0.17 x 10" 4 . See also Fig.® for uj. 



4 Conclusion 

In this paper the phase transition from quintessence to phantom era in 
quintom model has been discussed. By the assumption that the Hubble 
parameter has a Taylor expansion in terms of comoving time, the behavior 
of the Hubble parameter and the potential near the transition time, and the 
relation between them, have been studied. See eqs.([14 jl -([16 jl . The conditions 
need to satisfy to have a quintessence to phantom phase transition have 
been calculated, i.e. evenness of a and eq. (fl"8|) . The same results obtained by 
studying the relations between equation of state parameter uj, the quintom 
potential and matter density, see ea. (|28j) . 

To be specific, we have considered special cases by determining the quin- 
tom potential. For slowly varying exponential quintom potential, we have 
shown that the equations are consistent when a = 2 and hi > 0. In this way 
it has been proved that the phase transition occurs with a rate depending 
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Figure 3: H 2 /h^ as a function of r, obtained from eq.Q (with initial con- 
ditions mentioned in the text) (points), and (1 + Hit 2 ) 2 (line). 



Figure 4: uj versus r, using approximation (|56|) (points) and H = /io(l + 
H\T 2 ) (line) in eq. (|l()|) for power-law potential. 
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on the density of matter at transition time, see the discussion after eq.(|48|). 
The relation of matter and dark energy fields and their densities have been 
obtained at transition time. See eqs. (|43j) . (|4"5|) and discussion after eq. (fl9|) . 
For slowly varying power law potentials, similar results has been deduced. 

In these examples we restricted ourselves to slow roll quintom model near 
the transition time. For j m ~ 1, we expect at the transition time H = 0, 
the kinetic energy of quintom field be of the same order of magnitude as 
dark matter density, see eq.®. Using this conclusion in eq.Q results that 
the main part of the energy density is provided by the potential energy of 
scalar fields. Therefore the dark matter density is very small with respect to 
the density of dark energy, hence the coincidence problem still unsolved in 
these examples. 
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